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singularly perturbed boundary value 
problems using nonpolynomial splines 
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Abstract 


In this paper, a class of second-order singularly perturbed interior layer 
problems is examined. A nonpolynomial mixed spline is used to develop 
the tridiagonal scheme. The developed method is second as well as fourth- 
order accurate based on the parameters. Error analysis is also carried 
out. The method is shown to converge point-wise to the true solution with 
higher accuracy. Linear and nonlinear second-order singularly perturbed 
boundary value problems have been solved by the presented method. Five 
numerical illustrations are given to demonstrate the applicability of the 
proposed method. Absolute errors are given in tables, which show that our 
method is more efficient than previously existing methods. 
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1 Introduction 


Second-order singularly perturbed boundary-value problems (SPBVPs) have 
gained more importance for two reasons. Firstly, they occur in many areas of 
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science and engineering, like combustion, control theory, nuclear engineering, 
fluid mechanics, chemical reactor theory, and so on. Secondly, the occurrence 
of sharp boundary-layers as ¢, the coefficient of the highest order derivative, 
approaches zero, makes it difficult for standard numerical methods. The ap- 
proximate solution to boundary-value problems (BVPs), in which the highest 
derivative is multiplied with a small positive parameter, is described here. It 
is a well-known fact that the solution to the SPBVP exhibits a multiscale 
character. There are a number of methods available in the literature for 
solving SPBVPs [8, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 18]. Some other recent 
methods for solving nonlinear problems are given in [15, 16, 17]. 

We consider second-order SPBVPs of the form 


ez) (x) = G(x, z, 2'), r<a<s, (1) 
with 


2(r)=™, — 2(8) =n; (2) 


where 0 < € < 1 is a parameter multiplied by the highest derivative, 
G(a,z,2') are continuously differentiable functions on [r,s], and no and ny, 
are constants. 

Consider the following assumptions, for r < x < s and —co < z < 2’ < oo: 
(i) G is continuous, 

(ii) 9S and 8S exist and are continuous, 

(iii) 9¢ > 0 and | ge |< D, for an arbitrary constant D. 


The above assumptions from [11] confirm the existence and uniqueness of 
the solution to the BVP (1). 
Here, we develop a method using a nonpolynomial spline for solving the gen- 
eralized form of second-order SPBVPs. We use three-point finite difference 
approximations, which gives in place of give the accuracy of order four. When 
we implement the method, a tridiagonal system is obtained and solved using 
MATLAB. 
The paper is organized as follows: In Section 2, the derivation of the numeri- 
cal method is presented. The obtained scheme is applied to the second-order 
SPBVPs in Section 3. Error analysis of the method is discussed in Section 
4. Numerical results are given in Section 5 to prove the efficiency of the 
proposed scheme. Finally, concluding remarks are presented. 


2 Nonpolynomial mixed spline 


Let the interval [r,s] is divided into n equal parts as r = a < a1 < 4 < 
+++ << 2p, = 8, by introducing 


u=rt+ih, whereh=(s—r)/n, i=0,1,...,n. 
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Let 
Ri(x) = aje®*-*) + b;[sin(6(a — 2;)) + cos(d(a — 2;))] + ci, (3) 


be a function defined on [r,s], which becomes an ordinary quadratic spline as 
parameter 6 —+ 0 and 6 > 0. 


To evaluate a;,b;, and c;, we let 


(S; + Si41), O<i<n. (4) 


Nlre 


Ri(ei)=%, Riltin) =x, RP (e) = 


Using the above conditions, we get 


S41 — % 1 e? —1 2 
a; ¢ 292 [ ¢ | (S; + Sif), 
_ St1 — Xi 1 e? -—1 2 ; : 
b; ; agel—- Ih (Si + Sixt), 
Zait(2+OQ)z%;, 1 28-1 ‘ 
: ; Tae ji (S824), 


where ¢ = —2+sin@+cos@ +e? and 6 = Xh. 
Using the first derivative continuity condition, R*_,(2;) = R¥(a2;), k = 0,1, 
the following scheme is obtained: 


bzi—1 + 62, + eign = h?(b1Si_-1 + Si + V1 S41), 1<i<n-1, (5) 


where 
sin @ + e? + cos0 
o= — 
2+sin@ — e? — cos6 
é= ; 3 
~=1, 
i (2sin 6 — 1)e® + cos@ — sind 
La 462 ? 
sin de? — sind 
oa 
20 
e? — sin@ — cos6 
= : . 
40 


Remark: Our scheme reduces to [2] when 


(¢,€, 0, 1, &1, 1) _ (12, 1, 


Nl Re 
Ble 


1 
4’ 
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Error 


After expanding (5), using the Taylor series, we obtain the truncation error 
(TE) as follows: 


= (G+Et Wat (6+ What + WIPE — (te tHje 
sere ee or nore (AP 0 
pO EE (ELEM A 4 nee — (SPL 
h'| ~ ¥_( S Py -o(ns), 1<i<n-1 6) 


For various values of coefficients obtained in (5), we obtain the second- 
order method as well as the fourth-order method. For ¢; + €; + ¥, = 1 and 
$1 = V1, we obtain the second-order method. The TE for (¢, €, , $1, &1, U1) = 
(1, -2,1,1/6, 4/6, 1/6) is given as follows: 


1 
ti = hel) +O(h), 1<isn-1. 


For (¢,€,¥, 61, 1, U1) = (1, -2,1,1/12,10/12,1/12), we obtain a higher or- 
der method, that is, fourth-order. The TE is as follows: 


1 
i= ~The +O(h"), 1<i<n-1. 


3 Implementation of the scheme 


We consider the second-order SPBVP of the generalized form (1) as 
ez) (x) = G(a, z, 2’), r<ac<s. 
In particular, we take the linear second-order SPBVP as follows: 
ez) (x) = —p(x)2' (x) — q(x)z(x) + g(a), r<ac<s. (7) 


Now, implementing the scheme (5) on problem (7), we get the method, which 
is of second-order accuracy as follows: 


ebzy_1 + bz; + eben =h?(O1Gi1+4Git iG), 1<i<n-1, 
(8) 


where ¢; = Wy, and ¢, + €£; + ¥, = 1. The method of fourth-order accuracy 
is given as follows: 
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ehzi_1 + 62, + edzigy =h?(bGi1t+4GitwiGiu1), 1<i<n-1, 


(9) 
where (¢1, €1, 1) = (1/12, 10/12, 1/12) and 
G; = G(a, 25% Ar 
Gi_-1 — G(a, ety 2,4) 
Gi41 = = G(a, 241, Bs 
G,= G(a, 25% %), 


in which 


1 Zj41 — 4-1 


a ah i 
a — 32-1 + 42% — %41 
t—1 oh ’ 
af = 2-1 4z; + 3241 
i4+1 oh ’ 


s ai eg h 
Zz = cers : 20 (Gis4 Gi_1). 


Using the above approximations, we obtain an expression for the second-order 
method as follows: 


(ge — Shoipi-t = shea, - shia +h? dbigi—1) 21-1 

+ (€€ + 2hbipi—1 — Whar pisa + h7brgi)z% 

oh Ker shoapi-s ae shea: a Shubpiss + hei qin1)Zi41 
=W(dig-ithigtvigin), lsisn-1. (10) 


Similarly for the fourth-order method, it reads as follows: 


3 ih 1 1 Eih 
[eb 5 hor — Sy papi — Shérpi + Sha + 


Pi)Pi+4 


h 
+h? (b1 — ee Zi-1 


+ [e€ + 2hdipi—1 — 2h(y1 4 Ps 41 th brqi] x 
[ew shor an pi)pi-1 4 hep Sh(v Se piu 
+ Qn + 2 ot para }zes 
=h?[(o1 ae pi)gi—1 + gi + (Wr 4 a pian] l<ign-1. (11) 
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4 Convergence 


Here, we discuss the error analysis of higher order method, that is, fourth- 
order. In the matrix form, (11) can be written as 


MZ =U, (12) 
where M = (m,,;) is the tridiagonal matrix whose entries are given by (11) 
and U = [u1,u2,...,;Un—i]” is given as 
h h 
uy =h?[(o4 lp, )go + £191 + (¢ + p10 
3 h 1 1 
(ed 5dr lps )Po 5 iP + gh + wh )p2 
2 har 
+ h*(o1 —- 50? P1)40) Yo; 
h . 
uj =h?[(d1 — a ee + E19 + (4 Sty ocsal 2<1<n-2, 
20 20 
héy 
Un-1 =h?[(¢ — Pan 1)9n—2 + €19n—1 + (Yi + Mey 1)9n] 
1 3 
(ew 5 or MDa 1)Pa—2 + shea + ha + Sn 1)Pn 
hy 
+ h? (4 + mn Pn—- 1)9n)Yns 
and Z = [z1, z2,.--,2n—1]’. We also have 
MZ=U+4+T, (13) 
where Z= is Penn eel and T = [t1,to,...,¢ n- i]’. 


Also t; = [% — zi]?,i = 1,2,...,n —1, is the local TE. From (12) and (13), 
we have 


M(Z-Z)=T, 
ME =T, 
E=Z-4Z. 


Now we calculate the sum of elements of each row of the matrix M as 


Cane e€ + ep + 3h(o1 — BYpi)po + Shep 
MORIN ap yy a ti h2é,q, +h? i=1 
sh(dit Sy P1)P2 + h°Erg (Wi +" 50 h&1 pga, i=l, 


Ce ie ad 
/ +h? ((b1 — Bp )qi-1 + git (di + Spi) ai41), 2<i<n-2, 
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Sone get € 4 gh(da "30 Pr—1)Pn—2 — ShEPn—1 
—3h(bit+ 8 pn—1)pn +h? (1 — Bt pn—i)an—2 + h?Ergn—1, =n. 


Furthermore, 


3 1 : 
Sum; > h? (Eig + vi92 70 Si POP 40 S1P1P2) + O(h?), i=1, 
Sumi > h?(d1g-1+ fig + Yigigi), 2<i<n-2, 


3 1 : 
Sum; > h?(b1dn—2 + €1gn—1 — 7p otPuPn—1— Fo £1Pn—1Pn—2) + O(h3), i=n-1. 


Let 0 < K € Z* be the minimum of | p; | and | q; |. For sufficiently small A, 
we can say that 


Sum, > h?[(. + v1 — ik 
Sum; > h?[(¢d1 + & + v1) 4K], 2<i<n-2, 


Sumy. > U(r + & — EIR, 


Therefore, we get 


a 1 | 
Sum, ~ h?[(qoé1 + v1) K] 
u < : 2<I< 2 
’ StsSn— 4, 
Sum, ~ h?[(d1 + 1 + v1) K] 
il 1 
< 5 . 
Sumn—-1 ~ h?[(d1 + 761) 4] 
Furthermore, 
1 
1 Rete) KI)’ iI 
1 a 
Sum; $ ark Soe 
h?[(d1+ 4) K]’ ae 


The matrix M is monotone, irreducible and diagonally dominant for suiciently 
small h. Therefore, M~! exist and is positive. Hence, 


Zl] = ||e- "HIT. 


Let M~' = (mj ,). Then by [19], we get 
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n—-1 
~~ mj Sum; = 1, 1<j<n-1. 
i=1 : 
Therefore, 
o 1 
mM, * 
a Sum,’ 
n-1 n-1 1 
-1y = ane 
Inn l cape Imisl S 2, Sum; 
j= i= 
1 al 1 
= [ 2(9 a 2 2 9 ls 
(Ba tw)K Wdthaty)K h(gi+ pa)K 
n—1 
(= } <i<n-1. 
|Z; I oe \Til, l<i<n-l 
Finally, 


1 1 
patty bth ty dit wet 


Ell = |M~*I17 | s [ ale 


h2K 
As the fourth-order method, ||T'|| = O(h®) using (6). Then 


1 1 1 1 


E O(h®) = O(n’). 
| Sek faa eGo ae Oe. me) ve) 


Hence, the error is of order four. For the second-order method, ||T'\| = O(h*) 
using (6). Then || || = O(h?). We can prove the convergence of the second- 
order method by following the above procedure. 


5 Numerical illustrations and discussion 


To check the applicability of the developed method to existing problems, we 
solve three linear and one nonlinear problem of the type (1). The maximum 
absolute errors (MAEs) are tabulated in Tables 1-8. We have also solved the 
heat flow problem in Example 5. 


Example 1. Consider the nonlinear SPBVP from [18] as 
ez) (x) + 22' (x) +e?) =0, Osa <1, (14) 
with 
z(0) = z(1) =0. 


The analytical solution is 
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z(x) =In (2/0 + »)) eine, 


The results of our method along with comparison are given in Table 1. Also 
we have tabulated the MAE in Table 2. 


Table 1: Comparison of approximate solution of Example 1 with the method 
of [18] for « = 107% and n = 1000 


x Presented [18] Analytical solu- 
method tion 
0.001 0.5933418 0.6913641 0.5983404 
0.010 0.6835355 0.6825219 0.6831968 
0.020 0.6736751 0.6726859 0.6733446 
0.030 0.6639109 0.6629456 0.6635884 
0.040 0.6542413 0.6532992 0.6539264 
0.050 0.6446647 0.6437448 0.6443570 
0.100 0.5981091 0.5972949 0.5978370 
0.300 0.4309488 0.4304523 0.4307829 
0.500 0.2877780 0.2874905 0.2876821 
0.700 0.1625668 0.1624234 0.1625189 
0.900 0.0513068 0.0512663 0.0512933 
1.000 0.0000000 0.0000000 0.0000000 


Table 2: MAE of Example 


1 for ($1, €1, di) = (1/12, 10/12, 1/12) 
9-4 9-6 9-8 


n\e Q-2 
16 3.010 x 10-2 | 5.450 x 10-2 | 3.950 x 10-2 | 3.208 x 107! 
32 2.900 x 10-2 | 2.010 x 10-7 | 8.300 x 10-2 | 1.275 x 107? 
64 2.870 x 10-2 | 1.600 x 10-2 | 4.430 x 107-2 | 3.280 x 1072 
128 2.860 x 10-2 | 1.550 x 10-2 | 1.310 x 10-2 | 7.760 x 10-2 
256 2.860 x 10-2 | 1.540 x 10-7 | 5.300 x 10-2 | 4.150 x 10-2 
Example 2. Consider the linear SPBVP from [14, 18] as 
ez) (x) + 2!(2) = 1+ 22, O0<a<1, (15) 
with 
z(0) =0, z(1)=1 
The analytical solution is 
2 —1)\(1—e-*/e 
2(@) = a(@+1-— 2) 4 ia Le a 


(i —e-1) 


IJNAO, Vol. 12, No. 2, pp 483-497 


492 Khan and Khan 


The results of our method have been compared with the methods of [14, 18] 
in Table 3. We have also tabulated the MAE in Table 4. 


Table 3: Comparison of the approximate solution of Example 2 with the 


methods of [14, 18] for « = 10-3 and n = 1000 


x [18] [14] Presented Analytical 
method solution 

0.001 -1.0009970 -0.6311195 -0.6293169 -0.6298573 
0.010 -0.9918800 -0.9898546 -0.9878740 -0.9878747 
0.020 -0.9815600 -0.9796000 -0.9776400 -0.9776400 
0.030 -0.9710400 -0.9691000 -0.9671600 -0.9671600 
0.040 -0.9603199 -0.9584000 -0.9564800 -0.9564800 
0.050 -0.9493999 -0.9475000 -0.9456000 -0.9456000 
0.100 -0.8918000 -0.8900000 -0.8882000 -0.8882000 
0.300 -0.6114000 -0.6100000 -0.6086000 -0.6086000 
0.500 -0.2510000 -0.2500000 -0.2490000 -0.2490000 
0.700 0.1894000 0.1900000 0.1906000 0.1906000 

0.900 0.7098000 0.7099999 0.7102000 0.7102000 

1.000 1.0000000 1.0000000 1.0000000 1.0000000 


Example 3. Consider the two-parameter SPBVP from [20] as 


—eqz) (x) + €c2' (x) + 2(a) = cos(ra), 0<a<1l, (16) 
with 
2(0) = 0, x 1j=0 
The analytical solution is 
2(x) = Acos(rx) + Bsin(xx) + Ce™* + De~?4-*), 
where 
_ eqn? +1 _ Ee 
e224 (eqn? + 1)?’ 22 + (eqn? +1)?’ 
Table 4: MAE of Example 2 for (61, 1,1) = (1/12, 10/12, 1/12) 
n\e 2-2 Q-4 2-6 2-8 
16 8.8950 x 10-7 4.7391x10-* | 5.6800x10- | 4.688010 
32 5.5531 x 10-8 2.8359x10-° | 7.3000x10~-* | 2.2370 10— 
64 3.4736 x 1079 1.7529x10-© | 5.2469x10-4 | 5.8100 10~? 
128 2.1708 x 10-™ | 1.0925x10~” | 3.1398x10-° | 7.5000x10~ 
256 1.3640 x 1071! 6.8234 10-© | 1.9407x10-© | 5.3738x1074 
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Table 5: Comparison of the approximate solution of Example 3 with the 
method of [20] for n = 128 


Ec [20] Presented [20] Presented 
method method 
eq = 10-2 eq = 10-2 eq = 10-4 eq = 10-4 
10-3 9.2993x10-" | 2.6109x10~-® | 1.3294x10-* | 2.7811x10-4 
10-4 1.1557x10-* | 2.5996x10-® | 3.6708x10~4 | 2.7239x10~4 
10-5 3.4933x 10-8 | 2.5984x10-® | 2.8085x10-* | 2.7148x10-4 
10-6 2.6878 10-8 | 2.5983x1078 | 2.7232x10-4 | 2.7138x10~4 
10-7 2.6072 10-8 | 2.5982x10-8 | 2.7147«x10-7 | 2.7139x10-4 
l+e*% 1+e" 
C= A aa! D = A-—xy 
X _ €e — V2 + 4€q eo + Vee + 4eg 
aed 2€q : a 2€a , 


The results of our method along with comparison are given in Table 5. 


Example 4. Consider the SPBVP from [14, 18] as 


ez) (x) + 2'(x) — z(x) = 0, <a< 1, (17) 

with 
z(0) =1, 2(1) =1 
The analytical solution is 
qq) Coo 

where 

ae —1+V/1+ 4e _ -1-V1+4e 

2€ , 7 2€ 


The results of our method have been compared with methods [14, 18] in Table 
6. We have also tabulated the MAE in Table 7. 


Example 5. Consider the heat flow problem of viscous incompressible fluid 
over a stretching plate from [1] as 
2) (¢) + P.[1 — e~*]2! (2) + €[z(x)z (x) + [2'(2)]?] =0, 0 < a < 1,(18) 


with 
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Table 6: Comparison of approximate solution of Example 4 with the methods 
of [14, 18] for e = 10-3 and n = 1000 


x [18] [14] Presented Analytical 
method solution 
0.001 0.6003270 0.6007916 0.6011354 0.6007918 
0.010 0.3712379 0.3716054 0.3719728 0.3719724 
0.020 0.3749439 0.3753111 0.3756784 0.3756784 
0.030 0.3787160 0.3790830 0.3794502 0.3794502 
0.040 0.3825260 0.3828929 0.3832599 0.3832599 
0.050 0.3863742 0.3867410 0.3871079 0.3871079 
0.100 0.4062043 0.4065697 0.4069350 0.4069350 
0.300 0.4962382 0.4965853 0.4969323 0.4969323 
0.500 0.6062278 0.6065307 0.6068334 0.6068334 
0.700 0.7405963 0.7408182 0.7410400 0.7410400 
0.900 0.9047471 0.9048374 0.9049277 0.9049277 
1.000 1.0000000 1.0000000 1.0000000 1.0000000 


Table 7: MAE of Example 4 for (41, 1,71) = (1/12, 10/12, 1/12) 
n\e ie a a-8 g-8 
16 1.9745 x 10-© 4.1185x10-4 | 3.8100x10-2 | 2.9880 107! 
32 1.2290 x 10-7 2.4523x10-® | 5.0000x 10-2 | 1.4300x10-! 
64 7.6864 x 10-9 1.5138x10-° | 3.6087x10-7 | 3.7300 10-2 
128 4.8037 x 10-%° | 9.4320x10-® | 2.1567x10-® | 4.8000x 10-3 
256 3.0033 x 10-™ | 5.9003x10~-9 | 1.3326x10-5 | 3.4706 10-7 


where P, is the Prandtl number, which can take different values. Problem 
(18) in [1] is perturbed by the small parameter ¢ with the first part consisting 
of the terms independent of € and the second part involving ¢. Here, we take 
the first part of (18) as 

z)(¢) + P,[1 — e~*]z'(x) = 0, 


0<a<1, (19) 


with 


The MAE of (19) is calculated using the double mesh principle due to the un- 
availability of an exact solution. The MAEs are given in Table 8 for P,.=1.5. 
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Table 8: MAE of Example 5 


n 16 32 64 128 256 512 1024 
Second-order 2.1451) 2.7598x) 3.4967x| 4.3996) 5.5172») 6.9075x) 8.6421x 
method for | 107° 107-6 10-7 10-8 10-9 10-19 | 197-14 
(1, €1, Y1) = 

(1/6, 4/6, 1/6) 

Fourth-order 7.1156) 9.5771x| 1.2409x| 1.5788x| 1.9909x) 2.4996) 3.1305x 
method for | 107° 107-7 10-7 10-8 10-9 10-19 | 19714 
(1, €1, Y1) aa 

(1/12, 10/12, 1/12) 


6 Conclusion 


In this paper, we solved second-order SPBVPs using a tridiagonal scheme 
obtained by a nonpolynomial spline. The presented method is based on a 
lower degree spline that could reduce computational time and is more feasible 
working. The error obtained in the developed scheme (5) is of order two as 
well as four depending upon the choice of parameters. Point-wise convergence 
and MAE are given in Tables 1-8. Here, linear, nonlinear, and two-parameter 
perturbed problems are being solved. Comparison with the result of [14, 18, 
20] was given in Tables 1,3,5, and 6, which proves the efficiency of the method 
and also shows that our method is better than these existing methods. 
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